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By means of ab initio and tight-binding calculations it is shown that 
metastable nitrogen clusters Ng (boats) can form quasi-one- dimensional 
ensembles in which the nearest clusters Ng are bound to each other by 
covalent bonds. Those ensembles are characterized by rather high energy 
barriers (~ 0.3 eV) that prevent the fission of the ensembles into isolated 
Ng clusters and/or N2 molecules. 

A nitrogen molecule N2 is known to be the most stable form of nitrogen at standard 
conditions. Contrary to the overwhelming majority of other chemical elements, the 
cohesive energy per atom, Ecohin) = Etot{n) — Ei (where Etot{n) is the total energy 
of an n-atom cluster N„, Ei is the energy of an isolated nitrogen atom) decreases 
with n, so it is energetically favorable for clusters N„ with n > 2 to decompose into 
separated N2 molecules and thus to lower their energy. However, recently it was 
shown theoretically [|T], |^, ^ that several nitrogen clusters, e.g., N4 and Ng, can exist 
in metastable states which correspond to the local minima of the potential energy 
as a function of atomic coordinates. Although metastable clusters N„ with n > 4 
have not been synthesized yet, they attract much theoretical attention as potential 
high energy density materials (HEDM). The energy Eacdn) = Ecohi'^) — Ecohin) 
accumulated in the N„ {n > 4) systems approaches the value of Eacc ~ 3 eV/atom 
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In order to make use of N„ (n > 4) clusters as new efficient and environmentally 
safe energy sources, one should be able to create bulk metastable nitrogen structures 
composed of N„ {n > 4) clusters. On the one hand, in order such structures could 
exist, the nitrogen clusters should be connected with each other by covalent bonds. 
On the other hand, those bonds should not be "too strong" in order the clusters 
preserved their individuality and the energy accumulated in such structures was 
not substantially lower than the energy accumulated in an isolated cluster, i.e., the 
clusters should form ensembles. Besides, there should be a sufficiently high energy 
barrier preventing the fission of ensembles into isolated nitrogen clusters and/or 
nitrogen molecules. 

An extensive theoretical study of Mailhiot et al. |Q suggests that nitrogen can 
form bulk polymeric metastable structures at high pressures. However, the poly- 
meric structures cannot be considered as composed of nitrogen clusters, i.e., as 
ensembles of clusters. In this paper, we present the preliminary results of numerical 
simulations of bonding between nitrogen clusters and show that boats Ng can form 
metastable quasi-one- dimensional ensembles in which the nearest clusters Ng are 
bound to each other by covalent bonds. 

The version of the self-consistent-field theory with 3-21G, 6-31G*, and 6-311G** 
basis sets in HyperChem was used in the present work. First, we studied the ge- 
ometry and energetics of isolated clusters N4 (tetrahedron) and Ng (cubane and 
boat). Those clusters are shown in Figs. 1-3. The accumulated energies are 
^acc(4) = (1.98^2.21) eV/atom for the cluster N4, Eacc{S) = (2.44^2.59) eV/atom 
for the cubane Ng, and Eacc{8) = (1.53 ^ 1.66) eV/atom for the boat Ng, depending 
on the basis set used. Our results agree with previous works ||l], ||, If. 
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We have carefully analyzed a lot of different mutual arrangements of nitrogen 
clusters. In both cases of N4 and cubane Ng clusters, we have found no metastable 
configurations. However, in the case of boats Ng, we have revealed that two boats can 
form a "cluster molecule" (Ng)2. It is shown in Fig. 4. The energy accumulated in 
this molecule, Eacd^G) = 1-79 eV/atom (3-21G basis set), exceeds slightly the energy 
accumulated in an isolated boat Ng. We have demonstrated by direct calculations 
that the molecule (Ng)2 is mechanically stable, i.e., all its oscillation frequencies are 
real, giving evidence that the configuration shown in Fig. 4 corresponds to a local 
minimum in the total-energy surface as a function of atomic coordinates. 

In order to study the formation of large ensembles of boats Ng, we employed 
a version of TBMD (tight binding molecular dynamics) method with transferable 
interatomic potential 0, |[ used by us earlier for carbon clusters |^, ||, |10|, and 



modified for calculations on nitrogen systems. Although the TBMD method is not 
as rigorous as ab initio approaches, it allows one to study relatively large systems 
as well as to calculate directly the lifetime of the metastable state and the height 
of the energy barrier. For N4 and Ng clusters and (Ng)2 "boat dimers", the results 
of TBMD simulations agree qualitatively with ab initio calculations. We have also 
shown that rather long quasi-one- dimensional ensembles of boats can be formed 
through attaching of new boats Ng to the (Ng)2 "molecule". The ensembles (Ng)^ 
with m up to 6 have been shown to be metastable with accumulated energy about 
90 percents of the energy Eacd^) accumulated in an isolated boat. The heights of 
the energy barriers preventing the fission of metastable configurations exceed 0.3 
eV. 

Finally, we note a recent success in recovering the non-molecular phase of nitrogen 



at ambient pressure at temperatures below 100 K ||12|]. Structure and physical 
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properties of a new phase are not understood yet. Its characterization remains an 
important issue. Presumably, it is polymeric. This can be considered as a first step 
towards synthesis of bulk metastable nitrogen structures, including "cluster phases" 
studied in this paper. 
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Figure captions 

Fig.l. Tetrahedron N4. 
Fig.2. Cubane Ng. 
Fig.3. Boat Ng. 

Fig.4. "Cluster molecule" (N8)2 composed of two boats Ng. 
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